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ABSTRACT 

V1408 Aql (= 4U 1957+115) is a low-mass X-ray binary with an orbital period 
near 9.3 hr, whose compact star is a black hole candidate. The system shows 
a large-amplitude orbital photometric modulation at optical wavelengths. We 
have obtained new optical photometry of V1408 Aql from which we derive the 
orbital light curve and an improved orbital ephemeris. We show that the orbital 
light curve can be reproduced by a model in which the accretion disk around the 
compact star is thin, axisymmetric, and uneclipsed. The secondary star is heated 
by X-rays from the compact star and the accretion disk. The orbital modulation 
is produced entirely by the changing aspect of the irradiated secondary star with 
orbital phase. Because the system does not eclipse, the fits of the model light 
curves are insensitive to the orbital parameters, allowing a wide range of orbital 
inclinations and mass ratios. 

Subject headings: accretion, accretion disks, stars: neutron. X-rays: binaries - 
X-rays: individuals (4U 1957+115) 

1. Introduction 

Low-mass X-ray binaries (LMXBs) are close binary stars containing a black hole or 
neutron star surrounded by an accretion disk fueled by mass transferred from a low-mass 
companion. V1408 Aql is an LMXB with an orbital period of P = 9.329 + 0.011 h, derived by 



ThorstensenI ( 119871 ) from a roughly-sinusoidal modulation he found in its optical light curve. 
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Thorstensenl (119871 ) attributed the modulation to the varying aspect o f the secondary star , 
which is heated by flux from the primary star and its accretion disk. iHakala et al.l (Il999[ ) 
obtained two nights of UBVRI photometry, covering one orbital period in total, and also 
saw an orbital variation, but the light cur ve was more a s ymme tric than sinusoidal. Noting 
that V1408 Aql became bluer when fainter, iHakala et al.l (119991 ) interpreted the minimum in 
the light curve as a grazing eclipse of a cool, red rim on the accretion disk by the secondary 
star. In their model the ratio of the rim height to disk radius is large, H/ R ~ 0.2 — 0.5, and 
the orbital modulation outside eclipse is caused by variations in the height of the rim, which 
obscures more or less of the brighter inner disk. Because it invokes an eclipse, the model 
implies a high orbital inclination, the best fitting model having an inclination near 75°. 

4U 1957+115, the X-ray counterpart of V1408 Aql, is a persistent X-ray source, showing 
neither large transient increases in brightness nor faint quiescent states. Its X-ray flux does. 



however, vary on all time scales from 10 Hz 



of four) (IWijnands et al.ll2002uRussell et al. 



1-2 p ercent rms amplitude) to years (by a factor 



2010f). The orbital period has not been detected 



at X-ray wavelengths but iLevine fc Corbetl ( l2006l ) reported a "marginal detection" of a 
periodicity at 9.3175 ± 0.0005 h, perhaps related to the orbital period; and lNowak fc Wilms 
( 119991 ) repo rted a 117-day period in the X-ray light curve, although this period appears to be 
evanescent ( IWijnands et al.ll2002l ). The distance to 4U 1957+115 is unknown and its mean 
X-ray luminosity consequently uncertain; but for a distance of 7 kpc the mean unabsorbed 
luminosity is lo g(Ly) ?^ 36.5 (erg sec~^, 1.5-12 keV), while for a distance of 40 kpc it is 
log(Lx) ^ 38.1 jRussell et al.lboioh . 



The X-ray spectrum of 4U 1957+115 has been measured numerous times. It is unusually 
soft, placing 4U 1957+115 in the region of the X-ray color-color diagram occupied by black 
hole binaries, including black hole soft X-ray transients (SXTs) in their outburst states 



Schulz et al. 


1989; 


McClintock & Remillard 


2003) 



be the only known black hole binary in the Galaxy that has not been seen to cycle between 
low and outburst states. The nearly featureless spectrum can be successfully fltted with a 
variety of models that include some combination of three components: a soft component from 
a single- or multi-temperature black body, presumably coming from the accretion disk; a hard 
component that is either a power law, with or without a cutoff, or a thermal C omptonization 



spectrum; and a broad emission feature near 6.5 keV from Fe K-hne emission (lYaaoob et al. 



19931 : iSingh et al.lll994l : iRicci et al.lll995l : iNowak k Wilmslll999HWiinands et al.ll2002h . The 



spectrum is softer when 4U 1957+115 is fainter, harder when it is brighter. 



Most recently iNowak et al.l (l2008l ) flt hi gh-quality Chandra, X MM-Newton, and RXTE 
observations of 4U 1957+115 with diskbb ( iMitsuda et al.l Il984j ) and kerrbb models (see 
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Li et al. 


2005; 


Davis et al. 


2006) 



the XSPEC package, fixing the inchnation at 75°, and 
found that these black hole models fit the observed X-ray spectrum well. The fits re- 
quired high inner-disk temperatures, kTj„ = 1.2 — 1.8 keV, and a low normalization factor, 
{Rin/DioY cos i = 5 — 25, where i?^ is the inner radius of the disk in km, Diq is the distance 
in units of 10 kpc, and i is the orbital inclination. There are degeneracies among the derived 
values for the mass of the compact object, the mass accretion rate, and the distance, but the 
best fitting models ranged from a 3 M0 black hole with a spin of a* > 0.82 at a distance of 
10 kpc to a 16 Mq black hole with a* ^ 1 at a distance of 22 kpc. Nevertheless, the question 
remains: Is the compact object in 4U 1957+115 really a black hole? The evidence for a 
black hole is hardly conclusive. A so ft spectrum does not neces sarily require a bla c k hole 



- many SXTs contain neutron stars ( iTanaka &: Shibazakil Il996[ ). and Nowak et al. (^2008 ) 
did not explore models for disks around neutron stars. Furthermore, lYaqoob et al.l (11993 ) 
argued that the high values of kTin and low values of cos i deduced from fits to the 
X-ray spectral energy distribution are more typical of neutron star systems than black hole 
systems. 

In this paper we present new optical photometry of V1408 Aql. We refine the orbital 
period and derive the mean orbital light curve. There are no eclipses in the light curve. 
We achieve g ood fi ts to the orbital light curve with a model similar to that proposed by 
ThorstensenI (jl987l ) in which the orbital modulation is due to the varying aspect of the 
irradiated face of the secondary star. We believe this model is more physically realistic than 
models that invoke a tall disk rim to produce the orbital modulation, but it increases the 
range of possible orbital inclinations and mass ratios. 



Observations and Data Reduction 



We observed V1408 Aql from 2008 June 2 to June 7 UTC and again from August 2 to 
August 7 with the Argos prime- focus high-speed CCD photome ter on the 2.1-m Otto Struve 
telescope at McDonald Observatory ( iNather &: Mukadanj 120041 ). We did not obtain data on 
June 6 nor on August 3 and 4 because of clouds. The observations were made through a 
filter equivalent to a combined BVR filter transmitting the wavelength range 4130-7385 A 
and with an exposure time of 10 seconds for all frames. The data were reduced using a 
combination of IRAF routines and an IDL pipeline (R. Hynes, private communica tion) . 
Data reduction for the Argos photometer is discussed in detail in iHynes et al.l ( 120041 ). We 
measured relative fiuxes by dividing the fiux from V1408 Aql by the total fiux from two 
nearby comparison stars in the frame. Figure [1] shows a stacked image of three BVR images 
of the field of V1408 Aql for a 30-second effective exposure. V1408 Aql is indicated with 
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vertical hash marks and t he two comparison st ars are circled . The se are the comparison 
stars numbered 6 and 8 in iDoxsey et aLl (119771 ): iMargon et al.l ( 119781 ) . 



3. The Optical Light Curve and the Orbital Period 

To search for orbital variations we folded the light curve on a linear ephemeris with a 
period of 9.331 h (see below). Figure [2] shows the phase-folded data from each night binned 
into 350 equal-width phase bins. The error bars shown in the figure are based on photon 
counting noise and do not include other noise sources, such as accretion disk flickering. 
Observations with errors greater than ±0.001 in relative flux (~2.5%) due to poor seeing or 
clouds have been discarded. The light curve has a pronounced, roughly sinusoidal orbital 
modulation. There are also flickering variations with an amplitude of ~ 5% of the median 
flux over an interval of ~ 0.15 in phase or ~ 1.4 hours. There are no obvious eclipses. 

The mean flux from V1408 Aql changes on time scales of days to months. The system 
was fainter in June, brighter during the first two nights in August, and fading on August 6 
and 7. The top panel of Figure [2] shows the fainter light curve, the middle panel shows the 
brighter light curve, and the bottom panel shows the transition light curves. The color and 
symbol combinations indicate the different nights of observations. The dashed line is a sine 
curve fitted to the June data (top panel) and plotted in the middle and bottom panels for 
comparison. The brighter state has ~ 30% more flux than the fainter state but the peak-to- 
valley amplitude of the sinusoidal variations is nearly the same fraction of the median flux 
in both states. 

The orbi tal light curve of y i408 Aql in 2008 was nearly identical to the light curve 
measured by iThorstensenI (119871 ) in 1985. In both years the variations of the mean light 
curve were nearly sinusoidal with, as we will show in section 5.2, peak-to-peak amplitudes 
near 23%. While it is possible t hat th e mean orbital light curve was truly different in 1996 
when observed by lHakala et al.l (119991 ). they constructed the 1996 light curve from data that 
were obtained on just two nights and covered slightly less than one full orbit with no phase 
overlap. It is more likely that the light curve was distorted by flickering and night-to-night 
variations of the mean flux, and is not a good representation of the mean orbital light curve. 
Since the variations of color with orbital phase are also affected by these problems, it is 
premature to conclude t hat V1408 Aql becom es bluer at fainter orbital phases. On the other 
hand. Figures 2 and 3 in iHakala et al.l (119991 ) do show that any color changes are small and, 
in the BVR passbands, consistent with zero. 



Russell et al.l (120101 ) did not detect the orbital modulation at all in photometry of 
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V1408 Aql obtained between 2006 and 2009, and concluded that the system had changed. 
Their sparse data set (eg, 35 data points in V spread over 3.1 years) is not, however, well 
suited for detecting an orbital modulation that has been diluted by flickering and night-to- 
night variations. Also, their densest data set, from 2008, by chance covers only half the 
orbit with consequent loss of sensitivity to orbital modulations. Our data, taken within 
a few weeks of thei rs in 2008, shows that the orbital modulation was clearly present. The 
Russell et al.l (120101 ) data does show that the night-to-night variations can have a total range 
of at least 0.8 mag, somewhat higher than the 0.6 mag we observed. 

To refine the orbital period we first separately calculated the medians of the fainter state 
and brighter state light curves and then subtracted the medians to place the faint and bright 
light curves at zero median flux. We t hen perfornaed a Lomb-Sca rgle (L-S) period search, 
excluding the transition light curves (ILombl Il976t IScargld Il982[ ). T he L-S periodogram 
in Figure |3] shows a series of peaks near the orbital period found by iThorstensenI ( 1l987l ). 
with the two highest peaks at 9.331 ± 0.003 h and 9.393 ± 0.003 h. All these peaks have 
vanishingly small false alarm probabilities 10^^°), reflecting the large amplitude of the 
orbital variation, and all are two-month aliases of each other introduced by the two-month 
separation of our observing runs. Because Thorstensen's data were taken over a single 12-day 
interval and are not subject to one- or two-month aliases, we have used his period to choose 
among the aliases in our data. His period, 9.329 ± 0.011 h, agrees with the 9.331 ± 0.003 h 
peak in our L-S periodogram to within a standard deviation and differs from the other peaks 
in Figure [3] by six or more standard deviations. The orbital ephemeris is, then. 



T = 2454621.86(02) + 0.38877(13)E, 



(1) 



where T{E = 0) is the time of maximum flux. We searched for, but did not find, any periodi c 



modulations at other frequencies in the optical data, in agreement with iThorstensen 



(119871 ). 



4. The X-ray Light Curve and X-ray Periodicity Search 

Figure m shows the one-day-averaged X-ray light curve of 4U 1957+115 from the RXTE 
All Sky Monitor (ASM)Q during 2008 May to August. The times of our optical observations 
are also shown in the figure. The X-ray flux was fairly constant in June when V1408 Aql was 
faint at optical wavelengths. The X-ray flux was generally higher in late July and August, 
but on those dates in August when we observed V1408 Aql, the X-ray flux was not any 
higher than it was in June. 



^Data provided by the ASM/RXTE teams at MIT and at the RXTE SOF and GOF at NASA's GSFC. 
http://xte.mit.edu 
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We searched for periodicities in the X-ray hght curve from 1996 January to 2008 October 
in dwell-by-dwell samples using both a least-squares sine curve periodogram and a L-S peri- 
odogram. The L-S periodogra ms in Figures O and E] are fully-generalized L-S periodograms 
( iZechmeister fc Kursterll2009l ). This method is an improvement over the original L-S peri- 
odogram because it includes the weights of individual data points and allows for non-zero 
means. The periodograms have been normalized assuming the noise is white. Figures [5] and 
|6]also show least squares sine fit periodograms for comparison. 

Figure \5\ shows the periodograms of the X-ray light curve at periods near the orbital 
period. The open circle with error b a rs ma rks a frequency of 2.5726 day~^, corresponding to 
the orbital period from iThorstenseru (jl987l ): and the dashed line marks a frequency of 2.5722 
day^^, corresponding to the 9.331-hour orbital period from our data. The horizontal dashed 
line is drawn at the L-S power level for which the false alarm probability is 99%. The dot 
at 2.5 758 day~^ marks the "marginal detection" at P = 9.3175 hours from lLevine fc Corbet 



( I2OO6I ). There are no significant peaks in the periodogram at this or any other frequency 
between 2.54 and 2.6 day~^. Thus, there is no evidence for any modulation of the X-ray flux 
at periods near the orbital period. Figure E] shows the low frequency periodograms. The extra 
power at frequencies less than ~ 0.01 d~^ is real in the sense that the probability it could be 
produced by white noise is <C 10~^. No ne of the individual lo\v -freq uency peaks correspoii d 
to any of the long periods mentioned by lNowak &: Wilmsl (Il999[ ) and lWijnands et al.l ( 12002| ). 
including those at 58.5 d, 117 d, and 234 d. We therefore attribute this excess power to 
non-periodic or perhaps transient periodic variations of the X-ray flux. 



5. A Model for the Orbital Light Curve 



5.1. PreUminary Considerations 

We modeled t he orbital light curve of V1408 Aql using our light curve synthesis program!^ 
for X-ray binaries (iBayless et al.ll2010[ ). The program assumes that the orbital eccentricity is 
zero, the primary star is a point source surrounded by an accretion disk, and the secondary 
star fllls its Roche Lobe. The accretion disk can have a complicated geometry. It can be 
non-circular, non-axisymmetric, and spotted, and can have structures that extend vertically 
out of the plane of the d isk, such as an interior torus and a tall disk rim. The program uses 
stellar atmospheres from lKuruczl ( 1l996l ) for the spectrum of the secondary star and blackbody 
spectra for all other structures. It allows for heating by irradiation but not for second-order 



full description of the program is available at http : //pisces . as .utexas . edu/robinson/XRbinEiry .pdf . 
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heating by the irradiated surfaces. Light curves can be calculated for Johnson/Cousins 
filters or for a square bandpass over a specified wavelength range. Because we observed 
with a broadband BVR filter, we used the square bandpass option with wavelength range 
4130-7385 A. 



The X-ray spectrum of 4U 1957+115 is remarkably simple. iNowak et al.l (120081 ) suggest 
that it "...may be the cleanest disk spectrum with which to study modern disk models." 
As noted in the introduction, the continuum can be modeled with just a multi-temperature 
black body disk, perhaps with a weak optically-thin comptonizing corona. With the pos- 
sible exception of the Fe K-line, there are n o strong emission lines. The spectrum does 
have absorption lines, but iNowak et al.l (|2008[ ) attribute the lines to absorption in hot-phase 
interstellar material in the line of sight to 4U 1957+115. The X-ray light curve shows no 
periodicities that might be attributable to disk warping, super-humps, precession, or other 
disk asymmetries. Thus, neither the X-ray spectrum nor the X-ray light curve require mod- 
els with complex disk structures. We are also loathe to invoke a tall variable-height rim at 
the outer e dge of the disk to modulate the optical orbital light curve. As we have noted 
previously ( iBayless et al.l 120101 ). tall rims are highly unphysical, requiring either supersonic 
turbulence or temperatures much higher than expected near the edge of the disk. 

We have, therefore, adopted a thin, axisymmetric disk in our model. Since there are no 
eclipses, the only effect of the disk is to add constant flux to the light curve and to irradiate 
the secondary star; and the only source of orbital modulation is the varying aspect of the 
irradiated secondary star. This simple model has just five meaningful parameters: the mass 
ratio, the orbital inclination, the zero point in orbital phase, the X-ray flux from the central 
source, which heats the secondary star, and an additional constant optical flux from the disk, 
which also heats the secondary star. 



5.2. Fits to the Light Curve 

A preliminary analysis showed that the model yields synthetic light curves that fit the 
observed orbital light curves equally well for a wide range of parameters. We therefore 
abandoned any attempt to define "best" values for the mass ratio and orbital inclination 
and, instead, explored the range of parameters that yield acceptable fits. We considered 
three mass ratios, q = M2/MX = 0.3, 0.1, and 0.025, roughly corresponding to a compact 
object that is a neutron star, a low mass black hole, and a high mass black hole. For each 
mass ratio we fit synthetic light curves for models with inclinations of 10°, 20°, 50°, 65°, 
and 70°, plus some for models at 5°, 75°, and 80°. For concreteness, we set the mass of 
the secondary star to M2 = 0.4 Mq and its temperature to 3500 K. Our results are nearly 
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independent of M2 since the geometry of the system depends only on q and the scale of the 

1/3 

system depends only on . The secondary's low temperature prevents its un-irradiated 
surfaces from adding significant flux to the light curve. 



We quantified the quality of the fits with and used the simplex algorithm (INelder fc Mead 



19651 ) in IDL to find the values of the remaining three parameters (zero phase, X-ray flux, 
optical flux) that minimized x^- Because of the flickering, light curves from individual nights 
often depart from the mean orbital light curve by much more than the standard deviation of 
the measurements, which greatly increases x^- Since the flickering in effect adds non- white 
noise to the data and the light curves are sampled with a high cadence, the departures are 
strongly correlated, thwarting our attempts to artificially increase the variances of the indi- 
vidual data points to account for the flickering. As a result, the absolute value of loses its 
meaning. Nevertheless, the relative value of retains its utility for measuring the relative 
quality of fits. 

Figures [7] and [8] show the observed light curves overlayed by fitted synthetic light curves. 
The figures show fits to the light curves during both the fainter and brighter states, for two 
inclinations near the limits of the range we explored, i = 20° and i = 65°, and for the largest 
and smallest mass ratios we explored, q = 0.3 (Figure [7]) and q = 0.025 (Figure |8]). The 
peak-to-peak amplitude of the fitted synthetic light curve is 23% of the mean flux when 
V1408 Aql is in its fainter state an d 24% of the mean flux when it is in its brighter state. 



For comparison, iThorstensenI (119871 ) measured an amplitude of 0.116 mag or a peak-to-peak 
amplitude of 22.6%. The independence of the fractional amplitude from the mean flux is 
an important constraint. We interpret the constraint to mean that (1) the optical flux from 
both the disk and the secondary star is dominated by reprocessed X-ray flux and (2) the flux- 
weighted temperatures of the irradiated surfaces are high enough that the emitted optical 
flux is on the Rayleigh- Jeans tail of the blackbody distribution. Changes in the X-ray flux 
then produce proportional changes in the optical flux from the disk and from the heated face 
of the secondary star, leaving their relative contributions to the optical flux u nchanged. The 



blue colors of V1408 Aq l {B - V ^ 0.13, U - B ^ -0.70, E{B - V) ^ 0.3 flMargon et al. 



19781 : iHakala et al.lll999[ )) are at least roughly consistent with this interpretation. 



We achieved equally good fits to the orbital light curves for all three mass ratios. The 
dominant effect of increasing the mass ratio is to increase the relative size of the secondary 
star; but the effect of the increased size can be offset by reducing the amount of irradiative 
heating. As a result, the observed light curves do not by themselves yield useful limits on 
the mass ratio. 



We have two ways to place an upper limit on the orbital inclination. The first is from 
the lack of eclipses in the light curve. Let us assume that the accretion disk extends to its 
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tidal truncation radius, Rd = 0.9Rli, where Rli is the radius of a sphere with the same 
volume as the Roche lobe around the compact object. Eclipses of the disk would be detected 
in the mean light curve if they were more than about 2% deep. This occurs at i > 65° for 
q = 0.3, at z > 70° for q = 0.1, and at i > 75° for q = 0.025. The limiting inclinations are 
higher if R^ < 0.9Rli- The second way to place an upper limit on the inclination is from a 
flattening of the minimum of the synthetic light curves at high orbital inclinations. Figure [9] 
shows synthetic light curves for q = 0.3 (top, red line), q = 0.1 (middle, green line), and 
q = 0.025 (bottom, blue line), each for an inclination of 70°. The synthetic light curves all 
flatten at minimum light. This flattening is unrelated to eclipses but instead occurs because 
less than half the secondary star is irradiated, because the edges of the irradiated region 
are heated to lower temperatures, and because the secondary is not spherical. The lack of 
flat minima in the observed optical light curves places an upper limit of 70° on the orbital 
inclination. This second limit is independent of mass ratio and the size of the accretion disk, 
and is more reliable than the limits set by the lack of disk eclipses. 

The synthetic light curves for models with orbital inclinations of 20° and 65° both flt 
the observed light curves well. The flt for 20° has a slightly smaller but the difference 
is not statistically signiflcant. The dominant effect of decreasing the orbital inclination is 
to decrease the amplitude of the orbital modulation caused by the changing aspect of the 
secondary star, but the decrease of the relative amplitude can be offset by decreasing the 
fraction of the flux contributed by the accretion disk. This is shown in Table 1, which gives 
the fraction of the optical flux, a, that comes from the accretion disk as a function of orbital 
inclination; the fraction is also shown as the dashed lines in Figures [7] and [HI The value of a 
is nearly independent of the mass ratio and the luminosity state, but decreases with orbital 
inclination, dropping slowly from ~ 88% at i = 65° to ~ 67% at i = 20°, and then dropping 
rapidly to zero at ~ 5°. 

The flickering in the light curve disfavors inclinations less than 20°. The low values of 
a ai i < 20° would require an improbably high fraction of the disk flux to be modulated to 
produce the observed amplitude of the flickering. 

5.3. Constraints on q and i from Color and Temperature 

The results from the previous section show that neither the mass ratio nor the orbital 
inclination are strongly constrained by the shape of the light curve. They are, however, 
constrained by color and temperature, albeit rather weakly and with some uncertainty. We 
take as our starting point that the optical flux from V1408 Aql is dominated by regions of the 
secondary star and disk that are on the Rayleigh- Jeans tail of the black body distribution. 
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Let T2 and Td be the flux weighted temperatures of the irradiated areas of the secondary 
star and disk so that T2/Trf is the ratio of their emitted optical fluxes. We have two ways to 
calculate T^jT^- The first is from the amplitude of the orbital modulation, which is nearly 
independent of the luminosity and linear scale of the model; and the second is from the 
irradiative heating of the disk and secondary star, which depends strongly on the luminosity 
and scale of the model. The two calculations must agree for the model to be internally 
consistent. 

From the amplitude of the orbital modulation, the ratio of the optical energy coming 
from the secondary to the optical energy from the disk is 



Td J \Ad cos i J a 

where A2/2 and Ad cos i are the areas of the secondary star and disk that contribute signif- 
icantly to the optical flux, the factor of 2 because only half the secondary is irradiated and 
the factor of cosz because the disk is foreshortened. If the disk has expanded to its tidal 
truncation radius, then Rd = 0.9Rli, and equation [2] can be rewritten as 



T2 _ l-a 
Td a 



0.81(i?Li/a)2cosi 
0.5{RL2/ay 



(3) 



where a is the semi-major axis, and Rii/a and Rr,?Ja are funct ions of q and can be calculated 



from, for example, Eggleton's approximation ( lEggletonlll983l ). Table 1 gives the values of 
T2/Td implied by the fits to the light curves. Note that T2/Td must be larger for smaller 
mass ratios because the relative size of the secondary star is smaller. 

The temperature of the irradiated secondary depends primarily on q, the scale of the 
system (set by M2), the X-ray luminosity, and the fraction of the incident radiation that is 
locally thermalized and re-emitted. Table 2 gives typical values for T2 for X-ray luminosities 
of 10^^ and 10^^ erg s~^ for each of the three mass ratios, where we have assumed that the 
X-ray energy is produced near the compact object and that all the X-ray energy incident on 
the secondary star is thermalized and re-emitted. The temperature of an irradiated disk is 
extremel y sensi t ive to disk geometry and the location of the source of the irradiating X-rays 



(see, eg. 



HynesI (120051 )). As a result, our models give no reliable constraint on Td. We will 
simply take that Td must be greater than ~ 10^ K for its optical flux to be on the Rayleigh- 
Jeans tail of the blackbody distribution. It is possible, of course, that Td is substantially 
greater than 10^ K. 

We can now calculate T2/Td the second way, using T2 from Table 2 and taking Td > 10^ K 
so that the resulting T2/Td is an upper limit. For q = 0.3, the upper limits are T2/Td < 5 
for Lx = 10^*^ erg s~^ and T2/Td < 2.8 for 10^'' erg s~^ These values are consistent with the 
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values in Table 1 for any i > 20° and yield no additional useful constraint on the inclination. 
For q = 0.1 the upper limits are 4.0 and 2.2 and are mildly constrained, needing i > 30° for 
]^g38 gj,g g-i j > ^go ^q37 gj,g g-i_ At g = 0.025 the upper limits are 3.1 and 1.7 and 
become more strongly constraining. Only the model with Lx = 10^^ erg s~^ is internally 
consistent and even it requires i > 60°. 

These constraints are not, however, as certain as one might wish. The values of T2/Trf 
given in Table 1 are directly proportional to the adopted area of the accretion disk. If the 
radius of the accretion disk is as small as the circularization radius, T2/Td is reduced by a 
factor of 4 or more, and the inclination becomes essentially unconstrained even for q = 0.025. 
On the other hand, if Ta is greater than 10^ K, the constraints on the orbital inclination 
become more severe. Thus the most that can be concluded is that higher mass ratios are 
mildly favored and the models are internally consistent over a wide range of inclinations and 
mass ratios. 

Finally, it is possible that further observations m ay show that yi408 Aql does, indeed. 



becomes bluer at fainter orbital phases as claimed by iHakala et al.l (119991 ). Our model can 
reproduce such color variations - all that is required is that the heated face of the secondary 
star have a temperature less than the temperature of the accretion disk. This would strongly 
constrain the permitted mass ratios and inclinations, eliminating the q = 0.025 models 
altogether and restricting the g = 0.1 models to inclinations greater than ~ 65° and the 
q = 0.3 models to inclinations greater than ~ 40°. 



6. Summary and Discussion 



We obtained high-speed photometry of V1408 Aql on nine nights in June and August 
of 2008. The optical light curve flickered by at least 5% on time scales of minutes to hours 
and the orbit-averaged flux varied by 30% on time scales of a few days. The mean orbital 
light curve was roughly sinusoidal with a peak-to-peak amplitude n e ar 23% and was not 
measurably different from the light curve observed bv iThorstensenI (119871) in 1985. We 
attribute the lack of an orbital modulation r eported bvlRussell et al. and the different 

shape of the orbital light curve found by iHakala et al.l ( 1l999l ) to the sparseness of their 
data sets and to the distortions of the light curve induced by flickering and night-to-night 
variations. The reflned orbital period is 0.38877 ± 0.00013 d (equation [T]). There are no 
signiflcant periodicities in the RXTE ASM X-ray light curve at or near the orbital period. 



The orbital light curve is consistent with a model in which the accretion disk around 
the compact star is thin, axisymmetric, and uneclipsed. The secondary star is irradiated and 
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heated by X-rays from the compact star and the accretion disk. The orbital modulation is 
caused entirely by the varying aspect of the heated face of the secondary star with orbital 
phase. This simple model yields consistent fits for all mean light levels and it avoids an 
unphysically tall disk rim. 

At inclinations greater than 70° the synthetic orbital light curves flatten at minimum 
light. The lack of a flat minimum in the observed light curve places an upper limit of 
~ 70° on the orbital inclination. We cannot place an unimpeachable lower limit on the 
orbital inclination but argue that both the amplitude of the flickering and the lack of color 
variations over the orbit imply that the orbital inclination is greater than 20°. Thus, the 
orbital inclination lies between 20° and 70°. We achieved good fits to the light curves with 
mass ratios from 0.025 to 0.3. If the radius of the accretion disk is close to the tidal truncation 
radius, the smaller mass ratios, which correspond to systems in which the primary is a black 
hole, require X-ray luminosities nearer 10^® erg s~^ and orbital inclinations in the upper half 
of the permitted range. 



Figure 8 in iRussell et al.l (120101 ) plots optical luminosity against X-ray l uminosity for 



V1408 Aql and other LMXBs. The distance to V1408 Aql is unknown, but iRussell et al. 



(I2OIOI ) adopt three distances for comparison. If the distance is 7 kpc, the X-ray luminosity is 
~ 10^^'^ erg s~^, and at at 20 kpc, it is ~ lO'^^'^ erg s~^, placing V1408 Aql in a regions of the 
figure occupied by neutron star binaries but not black hole binaries. At 40 kpc, the X-ray 
luminosity is ~ iq^^-^ Q^g and V1408 Aql falls in a region occupied by some soft-state 
black hole binaries and the neutron star Z-sources. All these possibilities remain viable with 
our new model. 

If, however, V1408 Aql contains a black hole it would be unusual in two ways. First, 
as mentioned earlier it would be the only known black hole binary in the Galaxy that has 
not been seen to cycle between low and outburst states. Second, irradiative heating is more 
typical of neutron star X-ray binaries than black hole binaries. This is shown most clearly 
by SXTs during their outbursts. Large amplitude modulations from irradiative heating have 
been observed during outbursts of Aql X-1 and XTE 2123-058 and possibly X1608-52, all of 
which are classical SXTs containing neutron sta rs; and have also been observed in the related 



neutron star X-ray transient XTE J2 1 29+47 flThorstensen et all Il979l : IWelsh et all I2OOO 



Zurita et al.l l2000l : IWachter et al.l |2002| ). In contrast, orbital modulation from irradiative 
heating of the secon dary star is rare ly observed in black hole SXTs and, if it is observed, 
its amplitude is low (INeil et al.l 120071 ). This is due at least partly to the lower mass ratios 
of black hole binaries. At lower mass ratios the secondary star is comparatively smaller. As 
shown in section 5.3, the smaller secondary can produce a large-amplitude orbital modulation 
only if its irradiated face is heated to high temperatures. 
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Finally, the light curves of the black hole SXTs often show superhumps (iHaswell et al. 



2001; 


Zurita et al. 


2006; 


Kato et al. 


2009) 



binaries with disks whose radii cross the 3:1 o rbital resonance, which requires a mass ratio less 
than ~ 0.25 ( lLubowlll99ll ; iKato et al.ll2009l ). The failure of V1408 Aql to show superhumps 
in its light curve argues that its mass ratio is larger than ~ 0.25 and the mass of its primary 
star is less than ~ 1.6 Mp^. All t he known black holes in binary stars have masses greater 
than about 4 Mq (ICasaresll2006l ). While none of these arguments are conclusive, together 
they weigh heavily in favor of a neutron star and against a black hole in V1408 Aql. 



We thank R. 1. Hynes for helpful discussions. 
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Table 1: Properties of the Fitted Models 



Fraction of Flux from Disk Inferred^ 72/7d 



Inclination 


Faint State 


Bright State 


q = 0.3 


q = 0.1 


q = 0.025 


10° 


0.337 


0.334 


9.49 


24.9 


78.0 


20° 


0.683 


0.666 


2.28 


5.97 


18.7 


30° 


0.784 


0.759 


1.33 


3.48 


10.9 


40° 


0.830 


0.808 


0.88 


2.31 


7.23 


50° 


0.857 


0.839 


0.60 


1.56 


4.90 


65° 


0.880 


0.876 


0.29 


0.76 


2.38 


70° 


0.884 


0.869 


0.25 


0.65 


2.05 



^T2/Td is the temperature ratio required to produce the observed 
amplitude of the orbital modulation assuming — 0.9i?Li. 



Table 2: Typical Temperature of the Irradiated 

Surface of the Secondary Star 



Temperature (K) 
Mass Ratio Lx = 10^^ erg s"^ Lx = 10^^ erg s"^ 



0.3 
0.1 
0.025 



28,000 
22,000 
17,000 



50,000 
40,000 
31,000 
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Fig. 1. — A BVR image of V1408 Aql and its surrounding field. V1408 Aql is identified 
with vertical hash marks and the two comparison stars are circled. The image dimensions 
are 2^8 x 2f8 with north up and east to the left. 
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Fig. 2. — The binned light curves of V1408 AqL Each color and symbol combination indicate 
an individual night. Top Panel - The fainter state from 2008 June 2-7. Middle Panel - 
Brighter state from 2008 August 2 and 5. Bottom Panel - Transition from the brighter state 
to fainter state on 2008 August 6-7. The dashed hne is a sine curve fitted to the June data 
(top panel) and plotted in the middle and bottom panels for comparison. Because V1408 Aql 
was in transition on August 6 and 7, the orbital modulation is distorted, causing a spurious 
shift in the orbital phase of maximum light. 
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Fig. 3. — The Lomb-Scargle periodogram of the V1408 Aql optical data. The multiple 
peaks are all two-month aliases of each other introduced by the two-month separation of 
our observin g runs . The small circle and error bar mark the 9.329 ± 0.011 h period from 
ThorstensenI (119871 ). While the highest peak is at a f requency of 2.5552 d~^, the peak at 
2.5722 ± 0.0001 d"^ or 9.331 ± 0.003 h agrees with the iThorstensenI Jl987[ l period to within 
a standard deviation. It is this period that we adopt as the orbital period. 
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Fig. 4. — The one-day-averaged X-ray light curve from RXTE ASM. The times covered to 
the optical observations in Figure |2] are indicated with the vertical lines. 
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Fig. 5. — The weighted sine curve (top) and weighted Lomb-Scargle (bottom) periodograms 
of the RXTE/ASM X-ray hght curve of 4U 1957+ 1 15. T he open circle with error bars marks 
the orbital frequency determined by iThorstensenI (119871 ) and the vertical dashed line marks 
our refinement of the orbital period (equ ation [J) . The dot marks the "marginal detection" 
at P = 9.3175 h or 2.5758 d^^ found by iLevine fc Corbetl (120061 ). The horizontal dot-dash 
line in the bottom panel indicates the 99% false alarm probability. All of the peaks are 
consistent with noise. 
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Fig. 6. — The low- frequency portion of the weighted sine curve (top) and weighted Lomb- 
Scargle (bottom) periodograms of the RXTE/ASM X-ray hght curve of 4U 1957-1-115. 
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Fig. 7. — The observed orbital light curves of V1408 Aql in its low state (left panels) and high 

state (right panels). The observed light curves are overlayed with fitted synthetic light curves 
calculated for q = 0.3. The orbital inclinations are i = 20° (top panels) and i = 65° (bottom 
panels). The dashed line shows the contribution from the accretion disk, the remaining flux 
coming from the heated secondary star. 
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Fig. 9. — Synthetic light curves fitted to the observed hght curves of V1408 Aql in its fainter 
state. The orbital inclination is 70° for each light curve, and the mass ratios are q = 0.3 
(top, red line), q = 0.1 (middle, green line), and q = 0.025 (bottom, blue line). The models 
have been normalized to 1.0 at maximum and then offset by 0.25 for clarity. At this and 
higher orbital inclinations the minimum of the light curve is flat bottomed. 



